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1. INTRODUCTION 

Earth turbulent atmosphere introduces spatial and temporal variations in wave-front 
that lead to image degradation of optical systems. Astronomical telescopes, laser beam 
projection systems, and optical communication systems are limited by the presence 
of turbulence. In particular, the resolution of a ground-based telescope is notably 
modified. Generally, the telescope aperture is assumed to be smaller than the outer 
scale of the turbulence, so spatial frequencies of the turbulence with wavelength of the 
order of the aperture diameter, D, impart a random tilt on the incident wave- front. 
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This wave-front tilt translates to simple image motion at the image plane. It is the 
dominant atmospheric aberration across the telescope pupil. Statistical characteriza- 
tion of the image motion is of paramount importance because of its implications on 
the design of adaptive optics systems. 

In order to characterize temporally and spatially the statistics of the wave-front 
phase if several sensing methods have been used. 1 They use single, 2 double, 3 ' 4 and 
multiple 5 ' 6 (Shack-Hartmann) aperture sensors to measure the wave-front tilt. The 
centroid of the long-exposure images formed by each aperture is directly proportional 
to the slope of the wave- front across it. 

As usual, the phase structure function is 



where p', p G IR 2 , and (•) stands for the average using some unknown probability 
distribution. 7 Whenever a Kolmogorov developed turbulence is present, under the 
small perturbation and near-field approximations, the latter is turned into the widely 
known 8 



where r$ is the Fried parameter 9 linked to the spatial statistical properties of the 
refractive index, and C 2 is the phase structure constant, roughly near to 6.88. 

Interferometric measurements have corroborated the expression in Eq. (J2J). Many 
of these measurements have been made under the conditions mentioned above. But, 
significant departures from the 5/3 exponent have been experimentaly observed. 2 ' 5 ' 6 ' 10. 

In particular, for near to the ground measurements, exponents in the range 
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(1,5/3] have been determined experimentally. It is well-known that atmospheric tur- 
bulence is not always in its fully developed state, thus deviations from this simple 
model are likely — non-Kolmogorov turbulence. The phase structure function can then 
be generalized 6, 12 to include these results as follows, 



where (3 is the exponent associated with the phase spectrum, r 0jj a is the generalized 
Fried parameter and C 2 ^ is a constant maintaining consistency between the power 
spectrum and the structure function of phase fluctuations. If a Kolmogorov spectrum 
is chosen: (3 = 11/3, C 2 , ^ ~ 6.88 and r 0)/3 = r ; thus, we recover Eq. (J2J). 

In order to model turbulence-degraded wave-fronts by Kolmogorov turbulence, 
Schwartz et al. 13 have suggested that these are fractal surfaces described by a frac- 
tional Brownian motion (fBm) process with Hurst parameter 5/6 and a fractal di- 
mension equal to 13/6. Fractal properties are attributed to both the spatial and 
temporal behavior, and they are directly related through the Taylor hypothesis, or 
frozen turbulence approximation. The value of the Hurst parameter is in accord with 
the predictability of real stellar wave-front slopes. 14 Moreover, several algorithms for 
adaptive optics have been designed based on this statistical prediction. 13-15 In partic- 
ular, McGaughey & Aitken 16 have suggested that the source of predictability stems 
from the low-pass spatial filtering done by the wavefront sensor. These authors have 
modeled the temporal behavior of wavefront slopes by filtering the increments of a 
fBm with Hurst parameter equal to 5/6. 

On the other hand, several authors have simulated Kolmogorov phase fluctuations 
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over a finite aperture by using the random midpoint displacement algorithm. 17-19 This 
is also an algorithm for generating fractional Brownian motion surfaces — see Ref. [20, 
pp. 487-490 and 497-499]. Moreover, the fBm model was used as an alternative 
stochastic process to model the turbulent refractive index in lightwave propagation. 21 
In this paper the turbulent wave-front phase is modeled as a fBm. The result of 
Schwartz et al. is extended to include non-Kolmogorov turbulence. In Sec.|2]the gen- 
eral formalism for fBm processes is presented. This formalism is applied to determine 
the Strehl ratio in Sec. I3IAI and the angle-of-arrival variance in Sec. 13 IB I The latter 
result is compared against the well-known and widely used variance first obtained by 
Tatarski. 

2. FRACTIONAL BROWNIAN MOTION AND ITS ASSOCIATED NOISE 

Usually, natural phenomena behaving randomly are labeled as noises. These noises are 
characterized through the estimation of their power spectrum, ~W(v). 22-24 Empirically, 
an enourmously wide range of these spectra have been observed to follow power-laws 
proportional to \v\ _/3 , for some exponent (3. Better known as l/J^-type noises; 25 they 
are classified according to the value of the exponent, e. g. Ref. [23, ch. 3]. 

Since its first formalizations in the earlies 1900's (independently modeled by L. 
Bachelier and A. Einstein) the Brownian motion caught the attention of physicists. 
It is the most common representative of l// 2 -type noises; thenceforth, processes with 
such spectra are known as brown noises. On the other hand, the derivative of the 
Brownian motion is called white noise. The fact that it only can been defined as a 
distribution (in some probability space) is found in its tail-divergent power spectral 
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distribution, i. e. /3 — 0. Afterwards, any process between these two, with power 
exponent < (3 < 2, is referred as a pink noise. The last category is for those 
processes with 2 < (3 < 3, they are considered black noises. 

A (stochastic) process X(t) is self-similar with index H if, for any c > 0, 

X(t) = c H X(c~H); (4) 

that is, both processes are equal in distribution. The coloured noises are self- similar, 
with exponent H = (j3 — l)/2, as the generalized Fourier transform of their spectra 
can show. 26 Also, it suggests the presence of a slowing decaying auto-correlation, and 
thus of memory. 

Nevertheless, this colour classification is rather rough. It is insufficient knowl- 
edge of the power spectra to create stochastic processes modeling the randomness 
of quantities observed in the real world. As these random quantities tend to appear 
in dynamics equations other properties are needed to give them sense, e. g.: biman- 
ual rhythmic coordination differential equation, 22 Black-Scholes market equation, 27 
ray-optics equation, 28 etc.. 

Because of all the 'good' properties it endows, stationarity is desired. It is said a 
processes X is (wide sense) stationary if 

(X(t)} = (X(0)) and Cov(X(t + r),X(t)) = Cov(X(r), X(0)), 

for any r G R. 

Unfortunately, self-similar processes are non-stationary: Var X(t) oc t 2H ^ const.. 
Since we have lost ergodicity, the Wiener- Khinchin theorem fails. On the other hand, 
the existence of stationary increments does not contradicts the self-similar property. 
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A process with stationary increments is such that the probability law of its increments 
X{t + t) — X(t) is independent of t. 

Natural phenomena exhibit in general a non-gaussian behavior. Nevertheless, it 
is usual to append these to a gaussian distribution since, in this way, they become 
analytically tractable, Ref. [23, p. 35]. Moreover, choosing this distribution leaves un- 
affected the memory properties described by the spectrum. That path will be followed 
here. 

There is only one family of processes which are self-similar, with stationary incre- 
ments, and gaussian: the fractional Brownian motion (fBm). 29 The normalized family 
of these gaussian processes, B H , is the one with 30 

B H (0) = almost surely, E[B H (t)} = 0, 

and covariance 

E[B H (t)B H (s)] =\{\A 2H + \s\ 2H -\t- s\ 2H ) , (5) 

for s,t 6 I. Here E[-] refers to the average with gaussian probability density. The 
power exponent H is the Hurst parameter and its range is bounded. While the con- 
dition H > guarantees their (mean-square) continuity, H < 1 avoids degeneracy. 26 
Another more intuitive argument can be drawn. It is well-known these curves have 
fractal dimension equal to 2 — H. 31 Because they are embedded in the plane, H > 0. 
On the other hand, continuous parameterized curves should have dimension greater 
than one, and thus H < 1. 

These processes exhibit memory, as can be observed from Eq. (jSJ), for any Hurst 
paremeter but H = 1/2. In this case successive Brownian motion increments are 
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as likely to have the same sign as the opposite, and thus there is no correlation. 
Otherwise, it is the Brownian motion that splits the family of fBm processes in two. 
When H > 1/2 the correlations of successive increments decay hyperbolically, and 
this sub-family of processes have long-memory. Besides, consecutive increments tend 
to have the same sign, these processes are persistent. For H < 1/2, the correlations 
of the increments also decay but exponentially, and this sub-family presents short- 
memory. But since consecutive increments are more likely to have opposite signs, it 
is said that these are anti-persitent. 

Fractional Brownian motions are continuous but non-different iable processes (in 
the usual sense), and only give spectra exponents between 1 and 3. Nevertheless, fBm 
processes can be generalized to allow derivatives. A simple dimensional inspection 
suggests that the latter should have spectral exponent equal to (3 — 2H — 1; thus, 
covering the range — 1 < (3 < 1. 

Formally, continous processes are not called noises since they can be integrated 
pathwise. That is, given any continuous process X(t) the integral 



exists for any realization of the integrands. On the other hand, noises are not pathwise 
integrable. That is, the integral 



is not the limit of area approximating sums for any realization, i. e., there is no 
calculus in the classical sense. 
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The first construction of a Stochastic Calculus was made by Ito around 1940 
for Brownian motions. Later, these results were extended to more general processes — 
semi-martingales 32 and infinite dimensional Wiener (Brownian) processes. 33 The White 
Noise Analysis due to Hida, focused in the white noise rather than the Brownian mo- 
tion as a fundamental entity, is of particular interest here. 

Like common sense suggests the lack of conventional derivatives should be over- 
come through distributions. This is the basic idea underlying the white noise calculus. 
The problem is thus to embed these distributions into the right probability space. Let 
<p be an element of the Schwartz space S(M) (the space of rapidly decreasing smooth 
real- valued functions), and u is any element of the dual <S*(K). Therefore, the white 
noise is defined as the bilinear map W such that, W(<f>) = (0, u>) = (0, W(uj)) — where 
(•, •) is the bilinear map and (•, •) the usual internal product in L 2 (M). The space <S*(R) 
turns out to be a gaussian probability space and its elements uj the events. Moreover, 
the pairing coincides with the Ito integral, i. e. ((f), u) = f R <ftdB, u and using its 
properties (Ref. [34, p. 15]) it is found: 



that is, the white noise, as it was defined, is the derivative of the Brownian motion. 

In the last decade different approaches have been given to extend the stochastic 
calculus to fBm. The range of persistent processes have been particularly fruitful, 35-38 
not only because of its applications in practical problems also for its regularity prop- 
erties. Duncan et al. 39 successfully extended the white noise calculus to this range by 
means of a tool termed Wick product. These ideas were recently picked up by Elliott 
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and van der Hoek 27 who have given a complete calculus for all values of the Hurst 
parameter. A brief outlook based on their work is given in Appendix A. 

3. WAVE-FRONT MODELING AND APPLICATIONS 

Let ip be the phase difference between the average and perturbed wave- front. As it 
was argued at the introduction, it is a realization of a fractal surface. Moreover, the 
small perturbation and near-field approximations guarantee structures functions like 
Eqs. (0) or (jSJ). That is, it has stationary increments. Now, as always, it is assumed 
that the process ip is gaussian — see for example Ref. [8, p. 293]. Its power spectrum is 
also observed to follow a power law; thus, it is self-similar. At least it is valid within 
the inertial range, which is limited by two characteristic scales — the outer and inner 
scales, L and Iq respectively. 



be the isotropic fractional Brownian motion (ifBm). It is gaussian, self-similar and, 
under condition (jB.3|) (given at the Appendix B), has stationary increments. There- 
fore, we can define the generalized phase difference as 



where C v is defined as in Eq.(J2J), and H = 5/6 in the Kolmogorov turbulence case. 
Its structure function is, 



Now, let 



B H (p) :=B H (\\p\\) = B H (p) 



(6) 




(7) 
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where the last step is made under the condition |(p — p')/r \ 3 ^ 2 <C 1, which guar- 
antees this process has stationary increments — see Appendix B. Observe that since 
the Fried parameter can be interpreted as the diameter of the coherence area of the 
perturbed wave-front the last restriction is compatible. As it was stated earlier at the 
introduction, the structure function power exponent is restricted to the range (1, 5/3] 
for near to the ground measurements. Therefore, the Hurst exponent is confined to 
1/2 < H < 5/6. 



A. Strehl ratio 

The Strehl ratio, S, is used as a figure of merit in order to characterize the quality of 
beam propagation and the performance of a telescope. It is defined as the ratio of the 
central irradiance of the aberrated point spread function and the diffraction-limited 
point spread function. 40 

Considering a circular aperture of diameter D receiving an optical signal it was 
shown: 12 



S 



16 



7T 



u du 



(cosu) 1 — uyl 







exp 
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where D w (p) is the wave structure function. In the near- field approximation the wave 
structure function is replaced by the phase structure function. By using Eq. (jHJ) leads 
to 

if? r 1 r n r C 2 /n\ 2H 

(9) 
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In the case of Kolmogorov turbulence (H = 5/6) the following well-known expression 
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is recovered: 41 



S — — f udu 
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Remember that, for small phase aberration a 2 << 1, the Strehl ratio can be expressed 
as a function of the phase variance: 42 S ~ exp(— a 2 ). This formula implies that the 
normalized intensity is independent of the nature of the aberration and is smaller 
than the ideal unity value by an amount proportional to the phase variance. Under 
the definition in Eq. (J7|). it is 

^exp(-^p-). (10) 

Therefore, not only through the Fried parameter the quality of beam propagation can 
be set but also the Hurst parameter is relevant. 

B. Angle- of- arrival variance 

The path difference, or wavefront corrugation, of the wavefront surface from the av- 
erage plane is simply 

z(p) = ±<p(p), (11) 

as usual, A is the wavelength. Light rays are normal to the wavefront surface within 
the framework of Geometric Optics. The angle-of-arrival, at each normal plane, is 

27T OXi 

(i = 1,2). Following classical arguments, Ref. [8, pp. 334-336], the angle-of-arrival 
variance is given by 

<& = {%) + (%)= /^[^M+WfeH^A 2 /^ 2 ^^), (13) 
Jr Jr 
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where W^tV) is the power spectrum of <p(p). It should be stressed that in this ex- 
pression it is applied the Wiener- Khinchin theorem; thence, the phase is modeled as 
a stationary random variable. 

A divergent integral is obtained under the assumptions of Kolmogorov turbulence, 
small perturbation and near field approximation. In order to make it summable Rod- 
dier introduces a high and low frequency cut-offs, D~ l and Lq 1 respectively. It is a 
'more realistic expression' where the aperture diameter and the turbulence outer scale 
are involved. That is, 



We can obtain the following result, integrating the above equation and considering 
D < L : 

A more precise relation was given by Tatarski: 43 



where the proportionality coefficient is in radians squared units. However, coefficients 
ranging from 0.342, 44 ' 45 0.358, 3,4 to 0.365 46 have been given. It should be noted that 
these coefficients were obtained by using just only phase differences, so it is necessary 
that the wave-front remains unchanged over the whole aperture. Then, the pupil size 
must be smaller than the inner scale of the atmospheric turbulence, i. e. D < l . 
Now, according to Eqs. (JZJ) and (fTTj) . we have that 





(6.88/27r 2 )A 2 J D^/ 3 r - 5/3 , 



(14) 
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where C z = XC^jl-n . Therefore, the angle-of-arrival is 



9f(p) = ~q^ = -C z r^—^ Xi , (15) 

where W H is the fractional white noise as defined at Appendix A. The total variance 
of the angle-of-arrival is 

<h(p) = m + el] (p) = c^ H E[(w H ( P )f] . 

Let us calculate the fractional white noise variance using its chaos expansion, Eq. (jA.5|) . 
and the Wick product in Eq. (|A.7j) . Since, it is 

oo 

W H {p)W H {p') = W H (p) o W H (p') + MhUpWhUp'), 

n=l 

and considering statistical dependent variables are treated as if they were independent 
with respect to the average when Wick multiplied, E [W H o W H ] = E [W H ] -E [W H ] = 
• = — the last step is for the noise being a zero-mean gaussian variable. We have 
the following 

oo 

E[W H (p)W H (p')] = Y, M hUp)MhUp')- (16) 

n=l 

Finally, from Eq. (|15|) and the latter equation: 

oo 

<h(p) = c 2 z t^ h Y,m h Up)MhUp)- 

n=l 

The sum in the above equation can be analytically calculated. Consider 

where the Fourier transform 47 property of the Hermite functions was employed. There- 
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fore, 



oo 2 /• 

X2m h Up)MhUp')~ dvdv 



I I \l/2~H I /,1/2-H 



n=l 



expi(z/p + v p 



° 2 H f J J 'I 1 1/2— -ff I nl/2-H x , I /\ ■/ , / A 

= — / dvai 1 \u\ 1^1 d(z/ + z/)exp2(i/p + i/p i ) 

= — / di/ | 1 1_ " 2jF/ exp w(p — p), (17) 

for the intermediate steps one has to use the orthogonality and parity of the Hermite 
functions. First note from this equation that the covariance is stationary. But, if we 
set p = p' the angle-of-arrival variance 

2 ril -2H° 2 H f j I 1 1-2H 

° m ,H = C z r Q — / du \u\ 

is divergent! Then let us follow Roddier's idea and introduce an adequate cut-off to 
the above, 

^"(arp r ° n J dVV 

_ T(2H +l)smnH C% _ 2 
~ 2 2H n(l-H) 2vr 2 ' 

Observe that for H equal to 5/6, it is cr^ 5 , 6 = 0.452552 — where <r^ is the variance 
obtained by Tatarski. The cut-off d /a Roddier notably reduces the value estimated by 
Tatarski and others. As it was pointed out earlier, the scales considered throughout 
this paper are above the inner scale: in particular, D > Iq. Then, the difference 
between these variances is plausible. 

As we have seen, the removal of high frequencies is due to the finite size of the 
aperture. 48 In fact, since many scales are involved this filtering must be introduced 
in order to smooth out. Let us properly introduce this effect. 
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Define the smoothed fractional white noise as follows: given <p p (s) = <j)(s — p), 

oo 

Wf(p) = (<f> P , W H ) = ]T(0 P , Mh^HM (18) 

n=l 

is a noise built up by the contribution of each white noise with weighted function (f) p . 
Therefore, the variance of the smoothed noise is 

oo 



n=l 



Since (0 P , M H £ n ) = c H f dv . 



1/2-ff 



£ n exp — iz/p, it is 



oo „ 
n=l JR 



i/J| exp -iv{p - p ) 



'1 2 I |1-2H 



(P " P) 



where </>*(z/) 



acquires the form 



-v) have been used. Finally, the generalized angle-of-arrival variance 



n 2 _ r 2 fKZ r 2 r -2H g; 
a m,H - L/ z^ Z7TC H r & 



1 2 I \1-2H 



(0). 



(19) 



Observe that the function is a distribution-like function, it must satisfy the condi- 



tion: 



(f)(s) = V2^(p(0) = 1. 



(20) 



The natural election for <fi is the Fourier transform of a pupil with diameter D. 

Its normalized version according to Eq. (|2()jl is, 

"2 Ji(vD/2) 



7T (uD/2) 

Therefore, using Eq. (fT9~|) and the pupil filtering function, 
„* - r 2 r 2 r~ 2H ! dv \v\ l ~ 2H 2 J K vD l 2 ) 

_T(2H+l)T(H + l/2)T(l-H)smiTH C% _ 2H 



7r 3/2 2 2H-3 F ( i y + \)T{H + 2) 2tt 



(21) 
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Finally, it is c^j 5 / 6 = 1.04313 cr^. As we remove high frequencies the noise becomes 
more regular, and the wave-front variance approaches to that of Tatarski. 

4. CONCLUSIONS 

This paper introduces a stochastic process, the iffim, to model the turbulent wave- 
front phase. Not only it gives the right structure function for non-Kolmogorov turbu- 
lence, but also adds well-known statistic properties of the wave-front phase. Moreover, 
our model allows to extend results for two relevant optical quantities: the Strehl ratio 
and the angle-of-arrival variance. The expressions for these quantities depend on the 
Hurst parameter, thus on the dynamic state of the turbulence. 11 Remember that this 
parameter is related to the site location where the measurements are made. 

In particular, the expression obtained for the angle-of-arrival variance when H = 
5/6, Eq. (|21j) . is almost identical to the classical one found by Tatarski when high 
frequencies are filtered out. Nonetheless, for a Hurst parameter different from the 
one above a dependence with the wavelength appears. Up to now, it is unclear for 
us if such dependence exists for non-Kolmogorov turbulence. That is, if the Fried 
parameter is independent or not from the Hurst parameter. 

Using the formalism presented here a wider range of power spectra can be studied. 
Such as multifractal processes where the power exponent changes across frequency 
ranges. 

Also, asymmetric power spectra gives rise to self-affine surfaces; the phase tp is 
scaled differently depending on chosen axis. Thus, two Hurst parameters can control 
this behavior, and this formalism is applicable again. 
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Finally, since phase distortions of a wave-front transform into amplitude distor- 
tions in the wave cross sections, a similar analysis should be possible for the amplitude. 

APPENDIX A 

The purpose of this appendix is not to give a complete exposition of the calculus 
developed by Elliott and van der Hoek, 27 but an introduction of the tools used in this 
work. 

First, let Mh be an operator defined for any < H < 1 such that, 

hQ>(u) = c H \v\ l/2 - H ^(v), (A.l) 

where the hat stands for the Fourier transform, c 2 H = T(2H + I)sin7ri7, and the 
function is defined as in Sec. El The generalized fractional white noise is 

W H {(j>) = {M H <j>,u} = f {Mh4>) dB = [ (j)dB H - (A.2) 

Jr Jr 

thus, the extension to the Ito integral (with deterministic integrator) is again a bilinear 
map. Nevertheless, the test function <p must belong to L 2 H (K) = {<p : Mh4> £ -^ 2 (K)}. 
The fractional Brownian motion can be constructed as B H (t)(u) = (M^l( 0i t), uj)— 
where l( 0j &)(i) is the indicator function. 49 Two properties of the generalized noise are 
worth mentioning: E[(M#0, u)} = 0, and the isometry 

m[{M H ^u){M H (j>\u)] = (Mjj(j), Mjj(f)'). 

Afterwards, the stochastic processes subject to the same probability space are defined 
through what is called Chaos expansion. Shortly, any stochastic process X can be 
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written as the formal sums 



X{uj) =^c a H a \ 



to 



with c 2 a = K[XTi, a } ja\. Here it is defined a\ = aja^! • • • ot. n \, the factorial of the finite 
non-negative integer multi-index a. While 7i a {uo) = YYi=i H ai (.{£,i, to)) represents the 
stochastic component of the process, and it is build up through the Hermite functions: 

with H n the Hermite polynomials. These functions form an orthogonal basis satisfying: 

oo „ 

y2£n{x)£n{y) = 6(x - y), and / dx £ n (x)£ m (x) = 5 n>m . 
„=i Jk 

Particularly, consider that M H l[o,t](s) = Y^n=i( M H^[o,t], £n) Cn(s), then 

oo 

B H (t)(co) = (M H l m ,u) = ^ (M H l m ,Cn) (Zn,to). (A.4) 



71=1 



This is the chaos expansion for the fBm, since it is univocally defined and the (£ n , uo) = 
T~le„((£n,to)), please refer to Elliott and van der Hoek. Moreover, the operator M# is 
self-adjoint, so 

(M H l m ,Z n ) = (l [0) t ] ,M H Z n )= [ dsM H Z n (s); 

Jo 

the fractional white noise has the expansion 

, oo 

-B H (t) =^2M H UmnM = W H (t). (A.5) 

n=l 
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The stochastic integral in equation ()A.2|) is then rewritten as 



,i=l 



ds d>(s) 



n=l 



s)W H (s)ds. 



(A.6) 



The fractional Brownian noise is the integrated version of the generalized one. Also, 
the smoothed noise is defined: W¥{t) = f R 4>(s — t) W H (s) ds. 

Finally, it remains to be introduced the Wick product. Suppose X{uj) = J2 a a aHa(yj) 
and Y{u) = YlpbpHpif^) then the product is 



a+/3' 



Among the properties it has, maybe the most relevant is E[X oY] = K[X] ■ K[Y]. For 
the particular case X(u) = a + Y^Li a n'He n {^) and Y{uS) =b + Ylm=i a mHe m {^) is 



(X o Y)(u) = X(uj)Y(cu) - J2 a nK 



(A.7) 



n=l 



APPENDIX B 



The variance of the increments of an isotropic fractional Brownian motion is, usinj 
Eq. ©: 



E 



B H (p)-B«(p 



H / J- 



\P~P 



,2H 



(b.i; 



It will be shown that when p' and p are 'near' enough the latter turns into 



E 



B H {p)-B"{p> 



\P~P 



i\\2H 



(B.2) 
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It can be seen from Fig. IB. II that p — p' = 8g + 8 r and \p — p'\ 2 = 5 2 = \\S r \\ 2 . Thus, 
Eqs. (IB.lJl and (|B.2j) can be compared as follows 

|| p _ p '||2tf _ || (5e + 5r ||2H_ 5 2H 



5| + 25g ■ <5, 



Therefore, the condition 



5 2 

St 

5 2 



S 2 + 2S e ■ 5^ j 



- 1 



5 2 

must be fullfiled to allow the replacement made at Eq. (jB.2|) . Setting p' = p(cos 8, sin 9) 
and p = (p, 0), it is 8 2 e = 2p 2 (l — cos 6) and 28g ■ 8 r = 25p{ \ — cos 8). Supposing 9 <^ 1 
and 5 <C 1, 

^ 2 + 2<5 * • Sr - e 2 (C + ? V *V«r a , 



5 2 v^ 2 5, 

If it is set 5 ~ 5~ 2 5 2 , then the condition 

5 < 55 1/2 ~ 5 e < 1 (B.3) 

arises verifing Eq. ()B.2|) . 
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